Quantum-Mechanical Non-Perturbative Response 
of Driven Chaotic Mesoscopic Systems 
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Consider a time-dependent Hamiltonian Ti.{Q, P; x{t)) with periodic driving x{t) = Asm{Qt). It 
is assumed that the classical dynamics is chaotic, and that its power-spectrum extends over some 
frequency range < ui^i. Both classical and quantum-mechanical (QM) linear response theory 
(LRT) predict a relatively large response for Q < tJd, and a relatively small response otherwise, 
independently of the driving amplitude A. We define a non-perturbative regime in the {Q, A) space, 
where LRT fails, and demonstrate this failure numerically. For A > Ap^t, where Ap,t oc h, the 
system may have a relatively strong response for > Wci due to QM non-perturbative effect. The 
shape of the response function becomes A dependent. 



The wall formula for the calculation of friction in nu- 
clear physics 0, and the Drude formula for the calcu- 
lation of conductivity in mesoscopic physics, are just 
two special results of a much more general formulation 
of 'dissipation theory' ||^-|^. The general formulation 
of the 'dissipation' problem Q is as follows: Assume a 
time-dependent chaotic Hamiltonian HiQ, P;x{t)). For 
X = const the energy is constant of the motion. For 
non-zero V = x the energy distribution evolves, and the 
average energy increases with time. This effect is known 
as dissipation. Ohmic dissipation means that the rate of 
energy absorption ('heating') is d{H)/dt = fiV^, where 
H is defined as the dissipation coefficient. In case of pe- 
riodic driving x{t) = Asmiflt), one should replace V"^ 
by the mean square value i(v4r2)^, and the dissipation 
coefficient /x(0) becomes frequency dependent. For sim- 
plicity we assume that conservative work is not involved 
in changing x. 

In case of the wall formula, {Q, P) is a particle mov- 
ing inside a chaotic 'cavity', and x controls the deforma- 
tion of the boundary. Ohmic dissipation (in the sense 
defined above) implies a friction force which is propor- 
tional to the velocity, where /i is the 'friction coefficient', 
and fiV^ is the 'heating' rate. A mesoscopic realization 
of such system would be a quantum dot whose shape is 
controlled by electric gates. In case of the mesoscopic 
Drude formula, (Q, P) is a charged particle moving in- 
side a chaotic 'ring', and x is the magnetic flux through 
the hole in the ring. Ohmic dissipation implies Ohm law, 
where V = x is the electro-motive-force, /i is the conduc- 
tance, and fiV'^ is the 'heating' rate. For a mesoscopic 
realization of such system note that ring geometry is not 
important. One may consider a simple two dimensional 
quantum dot driven by a time-dependent homogeneous 
perpendicular magnetic field (For the latter geom- 
etry it is better not to use the term conductance while 
referring to the dissipation coefficient /i) . 

In the general analysis of the 'dissipation' problem one 
argues that due to the driving there is diffusion in en- 
ergy space. This diffusion process is biased because of 
its E dependence, leading to systematic increase of the 
average energy. This is the reason for having dissipation. 



Therefore we find convenient from now on to consider the 
diffusion coefhcient De as the object of our study. The 
relation between d{H)/dt and De constitutes a general- 
ization of the so called fluctuation-dissipation relation. 
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Fig.l. Upper diagram: The various V regimes in the theory of 
quantum dissipation for Unear driving x{t) = Vt. Lower diagram: 
The various {Q,A) regimes for periodic driving x{t) = As'm(Qt). 
Note the analogy with Fig. 5 of Ref. js] with x ^ A and V *-* Afl. 
The QM-adiabatic regime (including the regime A < Ac, but ex- 
cluding the narrow stripes of QM-resonances) is defined by having 
vanishing first-order probability to go to other levels. See the text 
for further explanations and definitions of Ac and Aprt . 

Ohmic dissipation is implied if oc or De oc A"^ 
in case of periodic driving. Such behavior can be es- 
tablished within the framework of classical mechanics 
|^,|| using general classical considerations 0|. The clas- 
sical formulation of the dissipation problem |^ can be 
regarded as a systematic scheme that justifies the use of 
classical 'linear response theory' (LRT). The precise con- 
ditions for the applicability of the classical LRT result are 
further discussed in |^,^, and we are going to mention 
later on what we call 'the trivial slowness condition'. We 
are interested in the quantum mechanical (QM) theory 
of dissipation. The traditional derivation j|] of QM LRT 
leads formally to the same result as in the classical anal- 
ysis fl . Therefore from now on we no longer distinguish 
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between the 'classical' LRT result and the 'quantal' LRT 
result and use just the term 'LRT result'. As a matter of 
terminology, it should be noted that the QM formulation 
of LRT, also known as Kubo-Greenwood formalism, is 
completely equivalent |^ to the well known Fermi golden 
rule (FGR) picture. 

So let us assume that the obvious classical conditions 
for the validity of the LRT result are satisfied^. Now the 
question is whether, upon quantization, there are addi- 
tional ?i-dependent conditions for the applicability of the 
LRT result^. In the traditional quantum mechanical lit- 
erature, as well as in the recent mesoscopic literature, the 
focus is on the consequences of having finite mean level 
spacing A. This leads to the identification of the QM- 
adiabatic regime (extremely slow driving), and to the dis- 
cussion of either the Landau-Zener mechanism |^ or else 
the Debye relaxation absorption mechanism |^ for dissi- 
pation, as well as to the discussion of QM-resonances. 
The main observation of is that there is another 
regime, the non-perturbative regime (see Fig.l), where 
QM LRT is not valid. As strange as it sounds, this does 
not imply a failure of the LRT result. On the contrary, 
another observation of is that the regime where the 
classical approximation applies, is well-contained in the 
non-perturbative regime, hence the LRT result becomes 
valid again because of quantal-classical correspondence 
(QCC) considerations. However, if the system does not 
have a good classical limit (as in RMT models) this 're- 
covery' of the LRT result is not guaranteed. Moreover, as 



we are going to discuss later, QCC consideration cannot 
exclude the possibility of having a relatively large quan- 
tal non-perturbative response whenever the LRT result 
is small in comparison. 

The outline of this letter is as follows: (1) We extend 
the theoretical considerations of to the case of peri- 
odic driving. (2) We give a specific example where the 
LRT result fails because of a quantal non-perturbative 
effect. (3) We comment on the issue of localization. (4) 
We discuss the role of QCC considerations in the the- 
ory. Based on the theoretical considerations, the reader 
should realize that the existence of the non-perturbative 
regime is not related to having finite mean level spac- 
ing A, but rather to having finite bandwidth Af, = huj^i, 
where uj^i is the dropoff frequency of the LRT response. 
In the context of mesoscopic physics this bandwidth is 
known as the Thouless energy. 

Given Ti.{Q, P; x) with x — const, we can define a 
fiuctuating quantity J^(t) = —dH/dx. The autocorre- 
lation function of !F{t) will be denoted by C(t). The 
power spectrum C{lu) is defined as its Fourier transform. 
The intensity of fluctuations is deflned as = (7(0), 
and it is convenient to define the correlation time as 
T^i = C(0)/C(0). We assume for simplicity of presenta- 
tion that the single time scale t^, completely characterizes 
the chaotic dynamics of the system: The power spectrum 
of the chaotic motion is assumed to be continuous, and it 
is non- vanishing up to the cutoff frequency uj^^ — 2'!t/t^i. 
We assume that C{lo) is vanishingly small for lj > oj^i. 
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Fig. 2. The response of a quantum mechanical system is displayed as a function of A and C The evolution is determined by the WBRM 
model Eq.(^). The units of energy and time and amplitude are chosen such that A = 0.5 and h = 1 and cr = 1 respectively. Left: Plots 
of Ds/A'^ versus H/oJci for few values of A. For small uj the plots coincide as expected from Eq.(^. As A becomes larger the deviations 
from Eq.(|l|) become more pronounced, and we get response also for f! > Wci- Right: Plots of De/Dq versus Ajy/h for few values of 
r2/a;ci- The LRT result Eq.{y) implies Ds/Dq = 1 for n/a;^ < 1 and De/Dq = for r2/a;ci > 1- The purpose of the horizontal scaling 
is to demonstrate that A^^t rather than Ac is responsible for the deviation from this LRT expectation. Each 'point' in the above plots 
is determined by a simulation that involves typically 35 realizations of the evolution until we start to see saturation due to dynamical 
localization effect. The typical time step is dt = 10"*. In each step we verify that the normalization is preserved to an accuracy of 0.01%. 
In order to eliminate finite size effects we have used a self-expanding algorithm. Namely, additional lOfe sites are added to each edge 
whenever the probability in the edge sites exceeds 10~^^. The diffusion coefficient is determined from the fitting SE(t)^ = const X . For 
sub-diffusive behavior (/3 < 0.86) we set De = 0. The preparation of each 'point' in the above plots requires ~ 4 CPU days on Alpha 
XPIOOO machine. 
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Consider the time dependent case x(t) = A sin (Hi). 
The LRT result for the diffusion in energy is 



(1) 



and we shall use the notation Dq = ji^{An)'^. The most 
transparent QM derivation of this result is based on the 
FGR picture. The energy levels of the systems are En, 
and the mean level spacing is A. The Heisenberg time 
is in = 27r?i/A. The transitions between levels are deter- 
mined by the coupling matrix elements {dTi./dx)nm- It 
is well known [ p^ that for reasonably small h this ma- 
trix is a banded matrix. The bandwidth is A;, = 2Trh/T^i, 
and the variance of the in-band elements is = i^/^h- 
It is common to define the QM system using the four 
parameters {A,b,a,fi) where b = A;,/ A. It is also use- 
ful to regard the semiclassical relations Td = 2Tifi/{h/S) 
and V = (27r?i/A)cr^ as definitions, whenever the classi- 
cal limit is not explicitly specified [as in Random Matrix 
Theory (RMT) models]. The FGR picture implies strong 
response if and only if hVl < Aj,, leading to 17 < uj^i as in 
the classical case. Using the FGR picture it is straightfor- 
ward to recover Dq — {7T/2){h/ A){aAfl)'^ in agreement 
with Eq.(Q). 

The trivial slowness condition for the applicability of 
classical LRT [|,| is Vt,, < Sx"^. Here Sx''^ is the para- 
metric change that leads to the breakdown of the lin- 
earization of 'H{Q,P;x + Sx) with respect to Sx. Upon 
quantization there are two other parametric scales that 
become important namely Sx"^"" and 5xp,^t. The for- 
mer is the parametric change which is required in order 
to mix neighboring levels, while the latter is the paramet- 
ric change required in order to mix all the levels within 
the band. Hence we define 



Ar = 6xr - ^ (X 
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The above parametric scales, of time-independent first- 
order perturbation theory (FOPT), manifest themselves 
also in the time-dependent analysis. FOPT gives the 
following result for the probability to make a transition 
from an initial level m to some other level n, 
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The total transition probability is p{t) — X]l^t("'k'^) 
where the prime imply omission of the n = m term. In 
the regime fl < one obtains 
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p{t) = -^„A^ 
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for < i < T,i 
for r,i < i < 
for 1/0 < t < 



while in the regime O > uj^i one obtains 
, , 1 ly A^ 



(1- cos(m))2 for < t < 

3/2 for T,, < t < t^ 



In both cases for t > we have recurrences, and there- 
fore p{t) < pita)- [One should be more careful near reso- 
nances: There tn should be replaced by 2Trh/S, where 6 is 
the detuning]. The necessary condition for applicability 
of FOPT at time t is that p(t') < 1 for any t' < t, which 
can be written as p([0,t]) ^ 1. The necessary condition 
for the applicability of the FGR picture is p{[0, t^]) ^ 1. 
This is the FGR condition [|| which guarantees the sepa- 
ration of time scales <C Tp,! . The FOPT breaktime r^^t 
is defined as the maximal t for which p([0, i]) < 1. Now 
we can define a non-perturbative regime by the require- 
ment p([0, r^]) > 1. It is straightforward to observe that 
the non-perturbative regime is contained in the region 



A > ^prt where 
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A 



< O < 
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(4) 



The location of the non-perturbative regime is illus- 
trated in Fig.l. For completeness of presentation we 
have also indicated the subregion in < where we 
have first-order response equal to zero. The condition 
is p([0, ooj) <g; 1 or equivalently p^ta) ^ 1. This region 
contains the QM adiabatic regime, including the region 
A < Ac, but excluding the narrow stripes of resonances. 

We have defined the location of the non-perturbative 
regime, but we did not yet give a suggestion how Eq.(|^) 
should be modified. Using RMT assumptions with re- 
gard to {dH/dx)nm, and inspired by related studies of 
wavepacket dynamics pM , we expect the result 



fu^) X Do 



(5) 



where Wpr = V/ {6x-p-,t/T^i), and C is a numerical constant. 
A detailed derivation of (||) will be presented in the fu- 
ture. (The crucial step is to argue that at Tp^t there is a 
crossover from ballistic behavior to diffusion in the sense 
of [^ ) . For periodic driving this result should be averaged 
over a period leading to De oc with a = 1/2. 

We wanted to give a numerical example that demon- 
strate the non-perturbative response effect. Evidently, 
the simplest is to consider a time-dependent version of 
Wigner's banded random matrix (WBRM) model, 



n 



En 



x{t) B 



(6) 



where Eq is an ordered diagonal matrix, and B is a 
banded matrix. This model ||l^,|l^ is characterized by 
the parameters (A, 6, a, K) which we have defined pre- 
viously. For the numerical experiment we have as- 
sumed rectangular band profile such that all the elements 
< |n — to| < b are taken from the same distribu- 
tion, and outside the band all the elements are identically 
zero. The results of the simulations are summarized in 
Fig.2. It should be reahzed that WBRM model Eq.(|) 
has a big disadvantage. Namely, unlike the physical ex- 
amples of the introduction, the statistical properties of 
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the model are not invariant for x{t) i-^ x{t) + const. 
One may wonder why we do not use one of the two 
other popular variations of Wigner model eg 
H = E + (cosx)Bi + (sina;)B2. The problem is that for 
these models one obtains dx^^t ~ Sx"^ ~ 2-k. Therefore 
there is no regime there where LRT fails because of quan- 
tal non-perturbative effect^. (Such failure requires the 
generic separation of scales feprt ^ ^x"^)- Thus it seems 
that the only way to make a RMT model x- invariant, is to 
keep it 'perturbative' in nature. The lack of a;-invariance 
in the standard WBRM model complicates the calcula- 
tion of the period-averaged and leads to cx A^~°' 
with a changing gradually from 1/2 to 1 . The numerical 
analysis of Fig. 2b fits well to a ~ 3/4. 

There are two types of localization effects that we had 
to consider in our numerical experiments. The "WBRM 
model localization" can be avoided by using amplitudes 
A < 6^/^(A/cr) in order to guarantee that the instanta- 
neous eigenstates of Eq.(^) are not localized at any time. 
The "dynamical localization effect" on the other hand 
cannot be avoided. It is associated with the periodic na- 
ture of the driving. Extending standard argumentation 
one observes that the eigenstates of the (one period) Flo- 
quet operator have localization length ^ x A, and that 
the associated breaktime is i* = ^ x (2^^/^l). The two 
must be related by 2Dj,t* = (CA)^ leading (in the LRT 
regime) to the result t* = 27^^ (A/ Ac)^ x t^. In all our 
numerical experiments the diffusion has been determined 
for times where dynamical localization is not yet appar- 
ent. Another possibility, which we have not used, is to 
add a small noisy component to the driving, such as to 
mimic the typical experimental situation of having de- 
phasing time much shorter than t* . 

It is not obvious that the non-perturbative behavior 
that is implied by RMT assumptions, and applies to 
RMT models, should apply also to Hamiltonians that 
possess a well defined classical limit. On the contrary, 
the same considerations as in |]l^ can be applied in or- 
der to argue that RMT considerations are not compatible 
with the QCC principle. Here we are going to explain the 
main idea, and to define our expectations. 

Taking h to be very small, it is obvious that eventually 
we shall find ourselves in the regime where A ^ . Let 
us consider the dynamics during a specified time interval 
< t' < <. The time t is chosen to be much larger than 
Tel. On the basis of QCC considerations, we should be 
able to make h sufficiently small such that the quantum 
evolution becomes similar to the classical evolution up to 
the time t. The classical analysis implies that during this 
time the stochastic behavior is established. Therefore 
having detailed QCC during the time t implies that the 
quantal Z?e can be approximated by the classical result. 
This leads to a contradiction with the RMT prediction 
Eq.(^ in the domain Q < a;^, but not in the domain 
n > Lv^i. We are going to further explain this last point. 

Denote the energy dispersion by SE^^{t), and the cor- 
responding classical result by 5E^i{t). For sufficiently 
small h it should be possible to make a leading order 



approximation 6E^^{t) « 5E^i{t) + TPg{t), with 7 > 0. 
In the < Wei regime the first term in this approxima- 
tion is dominant. On the other hand for > the first 
term gives a vanishingly small result for D^. Therefore, 
without any contradiction with QCC considerations, the 
second term becomes important. Therefore we may have 
in principle an enhanced quantal response for U, > uj^i. 

In conclusion, we have defined a non-perturbative 
regime in the (fi. A) plane, where LRT cannot be trusted. 
We have demonstrated an actual failure of LRT for a par- 
ticular (RMT) Hamiltonian. We believe that for generic 
chaotic systems the RMT mechanism for diffusion com- 
petes with the classical mechanism. The actual response 
of the system is expected to be determined by the pre- 
dominant mechanism. The study of this conjecture is the 
theme of our future studies. 
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